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Abstract 



We demonstrate that field theories involving explicit breaking of conti- 
nous symmetries, incorporate two generic classes of topological defects 
j^h. each of which is stable for a particular range of parameters. The first 

<— | \ class includes defects of the usual type where the symmetry gets re- 

stored in the core and vacuum energy gets trapped there. We show 
however that these defect solutions become unstable for certain ranges 
of parameters and decay not to the vacuum but to another type of sta- 
ble defect where the symmetry in not restored in the core. In the wall 
case, initially spherical, bubble-like configurations are simulated nu- 
merically and shown to evolve generically towards a planar collapse. 
In the string case, the decay of the symmetric core vortex resembles 
the decay of a semilocal string to a skyrmion with the important dif- 
ference that while the skyrmion is unstable and decays to the vacuum, 
the resulting non-symmetric vortex is topologically stable. 
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1 Introduction 

Topological defects]!], @ are stable field configurations (solitons ||) that arise 
during phase transitions [§J in field theories with spontaneously broken dis- 
crete or continous symmetries. Depending on the topology of the vacuum 
manifold M they are usually identified as domain walls § (kink solutions Q) 
when M = Z 2 , as strings and one-dimensional textures (ribbons [§, [7]]) 
when M = S 1 , as monopoles (gauged || ||, |10| or global 0, |12| ) and two di- 



mensional textures (0(3) solitons Ji"3"||3]])when M = S 2 and three dimensional 
textures |l4j (skyrmions||15||) when M = S 3 . Topological defects appearing 
in GUT phase transitions and quantum field fluctuations produced during an 
inflationary phase |1| |17|, [Tj| constistute two physically distinct mechanisms 



for the production of the primordial fluctuations that gave rise to galaxies 
and large scale structure in the universe. The macrophysical predictions of 
topological defects (particularly cosmic strings) for structure formation have 
been studied extensively [I], [19], |20|, [21], ^2], ^3], ^4| and compared with the 



corresponding observations with encouraging results. 

Defects corresponding to all of the above vacua are characterized by a 
conserved quantity called the 'topological charge' which gurantees their 
stability. It essentially counts the number of times the field configuration 
covers the vacuum manifold of the theory. This is achieved by the variation of 
the field at either spatial infinity or over the whole of space. Thus topological 
defects could be classified in two broad categories. In the first category 
the topological charge becomes non-trivial due to the behavior of the field 
configuration at spatial infinity. In this case continuity, forces the field to 
remain out of the vacuum manifold at a localized region in space where the 
spontaneously broken symmetry gets restored. This region is the core of the 
defect and is associated with vacuum trapped potential energy. Because the 
symmetric phase of the theory is retained in the core of these defects we will 
call them hereafter symmetric defects. Domain walls, strings and monopoles 
belong to this class. 

In the second category the vacuum manifold gets covered completely as 
the field varies over the whole of coordinate space. The field variable thus 
remains in the vacuum at all points of coordinate space. Its value at infinity 
is moreover identified with a single point of the vacuum manifold. Textures 
|fH|| (skyrmions [|15j), 0(3) solitons |TJ| (two dimensional textures J14], |25|j ) 



and ribbons || (one-dimensional textures [pqQ belong to this class which we 



will call for definitenss 'texture-like ' defects. Thus both the behavior of the 
field in the core of the defect (symmetric or non-symmetric phase) as well as 
its behavior at spatial infinity (covering or not the vacuum manifold) suffice 
to classify it as belonging to either the symmetric class or the texture class. 

Even though this classification can incorporate most of the well known 
defects it is not difficult to think of configurations that belong to neither 
class and yet possess quite interesting properties. These are configurations, 
appearing mainly in models with explicit symmetry breaking. There the 
field variable covers the whole vacuum manifold at infinity with the core 
remaining in the non-symmetric phase. Thus it is possible to have domain 
walls, strings and monopoles with cores in the non-symmetric phase of the 
theory. For definiteness we will call these 'non- symmetric ' defects. 

Examples of nonsymmetric defects have been discussed previously in the 
literature. Everett and Vilenkin [EBLin particular, pointed out the existence 



of domain walls and strings with non-symmetric cores. In their model how- 
ever the walls (strings) were bounded by strings (monopoles) formed at a 
phase transition at a higher energy scale. Thus, these defects were unsta- 
ble to shrinking and collapse due to their string tension. More recently, 



Dvali et.al. [27] and Bachas et.al. |7J considered a model which admits non- 
symmetric walls. They moreover showed the existence of classically stable 
bound states of such walls in this model which they called '27r walls'or mem- 
branes. They speculated that these states could be cosmologically disastrous 
even if the Z2 symmetry is anomalous. These bound state wall configurations 
were previously identified with one dimensional textures [Q or ribbons J| 
and are stabilized either by space compatification or by introducing a small 
anomalous term that breaks explicitly the Z2 symmetry [p8| , ||. 

A particular case of non-symmetric gauge defect was recently considered 



by Benson and Bucher [|29| (see also Refs. 0, pnfl ) who pointed out that 
the decay of an electroweak semilocal string leads to a gauged 'skyrmion' 
with non-symmetric core and topological charge at infinity. This skyrmion 
however, rapidly expands and decays to the vacuum. 

Further back and in the context of SU(5) GUT monopoles []TU[ higher 
SU(2) embeddings into the fundamental qunituplet space of SU(5) were 
shown to correspond to 'higher strength' monopole configurations |J2) . These 
were shown to be unstable and decay into the 'single-strength' fundamental 
monopoles while preserving the overall topological number. 

Our goal in this paper is to present more examples of topological defects 



that belong to what we defined as the 'non-symmetric' class. We will study 
in detail the properties of the simplest such configurations. More specifically 
we will first consider a model with a U(l) symmetry explicitly broken to Z 2 
which is inturn spontaneously broken. We will demonstrate that the model 
admits all three types of defects (symmetric domain walls, ribbons and non- 
symmetric walls) and identify semi-analytically the range of parameters in 
which each type of domain walls is topologically stable. 

In an attempt to study the cosmological effects of 2tt walls (ribbons) we 
perform numerical simulations of collapsing bubbles of symmetric and non- 
symmetric walls. We do not see evidence of the formation of bound states 
(27r walls). Instead, in the non-symmetric parameter range, the bubbles 
collapse asymmetrically (pancake collapse) and decay into the vacuum after 
a wall collision accelerated by the tension of the bubble. In the symmetric 
parameter range (symmetric phase in the wall core) the collapse is spherical 
and the bubble decays to the vacuum due to tension as expected. We give 
qualitative reasoning for this bubble behavior. 

In section 3 we extend our analysis to higher dimensional defects and 
in particular to the case of non-symmetric global strings. We identify the 
parameter ranges for stability of symmetric and non-symmetric strings in 
a theory with a global SU(2) symmetry explicitly broken to U(l). These 
results are then extended to the case of the semilocal string |33|, [54], ^, [36] 
where we identify the parameter range for the symmetric and non-symmetric 
phase. This result is in agreement with previous analyses |37], [36], |38j. We 



also briefly discuss the possibility of gauging the configurations considered 
which could lead to gauged strings and monopoles with non-symmetric core 



and thus to possible generalizations of electroweak strings |]39|, [40], |41], [|I| and 
monopoles. 

Finally, in section 4 we conclude, summarize and discuss the outlook of 
this work. 



2 Domain Walls: Symmetric vs NonSymmet- 
ric Core 

Consider a model with a U(l) symmetry explicitly broken to a Z 2 . This 
breaking can be realized by the Lagrangian density |l], p7| 

£ = l -d^d^ + -|$| 2 + ^-i?e($ 2 ) - ^|$| 4 (1) 

where $ = $1 + 2$2 is a complex scalar field. After a rescaling 

$ -> -^L$ (2) 

Vh 

x — >■ — x (3) 

m 

M ^ am (4) 

the potential takes the form 

V{*) = ~{o?\*\* + Re{&)-\\*\*) (5) 

For a < 1 this is a saddle point potential i.e. at $ = there is a local 
minimum in the $2 direction but a local maximum in the $1 (Fig 1). For 
a > 1 the local minimum in the $2 direction becomes a local maximum but 
the vacuum manifold remains disconnected, and the Z2 symmetry remains. 
This type of potential may be called a 'Napoleon hat' potential in analogy 
to the Mexican hat potential that is obtained in the limit a — > 00 and 
corresponds to the restoration of the S 1 vacuum manifold. 
The corresponding equation of motion for the field $ is 

$_ v 2 $ - (a 2 $ + $*) - |$| 2 $ = (6) 

which accepts the static kink solution 

$! = $ Jl = ±(a 2 + l) 1 / 2 taiih((^^) 1 / a x) (7) 

$2 = (8) 

This solution corresponds to a symmetric domain wall since the core of the 
soliton is symmetric ($(0) = 0) and the topological charge is resulting from 



(1.1 




V(<b) a<\/ 



Figure 1: (a) The domain wall potential has a local maximum at $ = 
in the $i direction, (b) For a > 1 (a < 1) this point is a local maximum 
(minimum) in the $2 direction. 



the behavior of the field at infinity (Q = ±($(-00) - $(+oo))/(o 2 + 1) 1/2 ). 
The form of the potential however implies that the symmetric wall solution 
may not be stable for a > 1 since in that case the potential energy favors 
a solution with $ 2 7^ 0. However, the answer is not obvious because for 
o > 1, $2 7^ would save the wall potential energy but would cost additional 
gradient energy as $2 varies from a constant value at x = to at infinity. 
A stability analysis is therefore needed and may be performed as follows. 
Consider an ansatz of the form 



$ = $ R (x) + 5§ 1 (x)e iuJlt + i8^ 2 {x)e lW2t (9) 

Substituting (9) in (6) we obtain to first order in the perturbations 



x 



6&1 + Stanh^-p)^!-^! 



cut 



s/2' 



or 



-<J$i 



x 



-6&2 + tanh 2 (^)<5$ 2 - 



o 2 -l 
o 2 + l' 



5$ s 



^2 

o 2 + l' 



-5$ 2 



(10) 

(11) 



For stability we demand that uof > i.e that there are no negative eigenvalues 
to the Schroedinger-like equations (10) and (11). The same conditions are 
obtained by perturbing the static energy functional 

irr,^ POO _ 1 

E = ^- dx(^' 2 + i-\^ - (a 2 \^\ 2 + Re(^ 2 )) (12) 

All JO I 

to second order around the solution (7), (8) and demanding that the pertur- 
bations contribute no negative part to the energy. 
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Figure 2: Field configuration for a symmetric wall with a = 1.65. 



The potential of the Schroedinger-like equation (11) is positive definite for 
a < 1 and it becomes negative and deeper for a > 1 as a increases. To find 
the critical value of a such that for a > a cr u (11) has negative eigenvalues 
we solve (11) numerically using the shooting method in the Mathematica E31 
system. The boundary conditions at the origin are <5$2(0) = 1 and 5<&' 2 (0) = 
0, corresponding to the ground state solution. We solve (11) for various values 
of c = ^Ti with uj 2 = looking for c crit such that lim^oo 5$ 2 (x) = which 
corresponds to a ground state with eigenvalue. Clearly for c > c crit the 
potential is deeper and there are negative eigenvalues. Using this method we 
find c cr u = §• This solution corresponds to a cr u = VS ~ 1.73. As expected 
a crit > 1 due to the effects of the gradient energy discussed above. Using the 
same method we also solved equation (10) and found that it has no negative 
eigenvalues as expected due to the topological stability of the component $ 1# 

We have verified this result using two other methods in order to both test 
the result and also test the accuracy of the methods. The first method is 
based on solving the full non-linear static field equations obtained from (6) 
with boundary conditions 



$i(0) 

$' 2 (0) 








lim $i(x) 



(a' + l 



|V2 



lim $o(x) = 



(13) 
(14) 



Using a relaxation method based on collocation at gaussian points 47 



to solve the system (6) of second order non-linear equations we find that 
for a < v^3 the solution relaxes to the expected form of (7) for $! while 
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Wall with Non — Symmetric Core, ex = 1.8 
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Figure 3: Field configuration for a non-symmetric wall with a = 1.8. 

$ 2 = (Fig. 2). For a > y/3 we find $i^0 and $ 2 7^ (Fig. 3) obeying 
the boundary conditions (13), (14) and giving the explicit solution for the 
non-symmetric domain wall. In both cases we also plot the analytic solution 
(7) stable only for a < \/3 for comparison (bold dashed line). As expected 
the numerical and analytic solutions are identical for a < a/3 (Fig. 2). 

The second method is based on a numerical minimization of the energy 
functional (12) using the steepest descent method of Refs. |7|, [0], |45|. We 
write the energy functional (12) on a one dimensional lattice of N points and 
a sum over 2(N — 1) variables 



E = ^ ^(i + i)-^) ^ + ( $ 2 (» + i)-$ 2 (») )2 + 



dx 



dx 



i($!(o a + $ 2 (*) 2 ) 2 - « 2 ($!« 2 + $ 2 « 2 ) - ($!(z) 2 - ^m 

We then evaluate numerically the quantities 

dE 



d^ii) 



and shift $,-(i) to 



where 



e ~ 10 -3 



(j = l,2), (i = l,...N) 



(15) 

(16) 
(17) 
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thus finding a new field configuration with lower energy. We repeat this 
procedure, keeping the boundary conditions, until the energy variation is 
negligible implying that we have reached a local minimum in configuration 
space. Our initial configuration was the analytic solution (7) perturbed by 
a small amount $ 2 = 0.01 e~ x . For a < a/3 the imposed perturbation 
decreased and the field configuration relaxed to the analytic solution (7). 
For a > v3 the perturbation increased towards the solution obtained by the 
relaxation method (Fig. 3). The verification of the result that a cr = v3 not 
only makes it strongly established but is also a test for the validity of the 
methods used. 

Bound states of pairs of non-symmetric walls have been discussed previ- 
ously in the literature under different names (one-dimensional textures [EH]. 



ribbons 0, 2n walls p7|). By Derrick's theorem |46| such configurations 



are normally unstable towards expansion induced by the gradient energy 



term \2~E\ but can be stabilized by introducing a small explicit breaking of 



the Z 2 symmetry [§, 27j or by imposing periodic boundary conditions 



The cosmo logical effects of these states were discussed by Dvali et. al. |P? 
They pointed out that the stability of a 2it wall system in the presence of 
an anomalous Z 2 symmetry can lead to a cosmic overabundance of walls and 
therefore restrictions can be imposed on the allowed range of parameters. As 
the cosmic horizon expands however, the 2ir wall systems will emerge as bub- 
bles of non-symmetric walls. It is the evolution and stability of such bubbles 
that needs to be studied in order to determine if these objects will cause a 
cosmological problem. 

In order to address this issue we constructed a two dimensional simulation 
of the field evolution of domain wall bubbles with both symmetric and non- 
symmetric core. In particular we solved the non-static field equation (6) 
using a leapfrog algorithm E7[ with reflective boundary conditions. We used 



an 80 x 80 lattice and in all runs we retained J| ~ | thus satisfying the 
Cauchy stability criterion for the timestep dt and the lattice spacing dx. The 
initial conditions were those corresponding to a spherically symmetric bubble 
with initial field ansatz 

$(t J ) = (a 2 + l) 1 / 2 tanh[(^^) 1 / 2 (p- Po )] +? 0.1 e Hl ^ Po1 ^- (18) 

where p = x 2 + y 2 and po is the initial radius of the bubble. Energy was 
conserved to within 2% in all runs. For a in the region of symmetric core 



Initial Non— Symmetric Wall Bubble, a = 3.5 
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Figure 4: Initial field configuration for a non-symmetric spherical bubble wall 
with a = 3.5. 

stability the imaginary initial fluctuation of the field $(£«) decreased and the 
bubble collapsed due to tension in a spherically symmetric way as expected. 

For a in the region of non-symmetric core stability the evolution of the 
bubble was quite different. The initial imaginary perturbation increased but 
even though dynamics favored the increase of the perturbation, topology 
forced the ImQ(t) to stay at zero on two points along the bubble: the in- 
tersections of the bubble wall with the y axis (Figs. 4, 5). Thus in the 
region of these points, surface energy (tension) of the bubble wall remained 
larger than the energy on other points of the bubble. The result was a non- 
spherical collapse with the x-direction of the bubble collapsing first (Fig. 
5). This asymmetric collapse of the bubble may be understood as follows: 
The bubble wall in a two dimensional projection may be approximated by 
a rectangle with dimensions x-y Let the x sides of the rectangle be in the 
symmetric phase and therefore have energy per unit length Eq. Then the y 
sides would be in the non-symmetric phase with energy per unit length E\. 

Clearly E > E\ since we are in the parameter range where the non- 
symmetric wall is stable. Thus the energy of the bubble walls is 

E = 2(E x + E iy ) (19) 



Evolved Non-Symmetric Wall Bubble, a = 3.5 
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Figure 5: Evolved field configuration (t = 14.25, 90 timesteps) for a non- 
symmetric initially spherical bubble wall with a = 3.5. 

which implies that it is energetically favorable for the bubbles to collapse in 
the x direction rather than the y direction. This argument clearly does not 
hold in the range of symmetric wall stability since in that case the energy 
would be E = 2E Q (x+y). Also it does not hold for an imaginary perturbation 
without the A factor since in that case the energy would be E = 2E\(x + y) . 
Generically, the two types of perturbations are equally probable and there- 
fore, in the parameter range of non-symmetric wall stability we anticipate 
about half of the wall bubbles to collapse asymmetrically (pancake collapse) 
while the rest would collapse in the usual spherical way. 

In all the cases we simulated we saw no evidence for formation of a bound 
state and in all cases the bubbles collapsed to the vacuum accelerated by their 
tension. This decay would be even faster for an anomalous Z 2 symmetry. 



Thus the cosmological problem discussed in Ref. [0 appears not to be 
realized in the cases we examined. 
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Vortex: with. Symmetric Core, ex = 2.6 
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Figure 6: Field configuration for a symmetric- core global string with a = 2.6. 

3 Generalizations to Other Defects 

It is straightforward to generalize the analysis of the previous section to 
higher dimensional defects like strings or monopoles. Consider for example 
a model with an £77(2) symmetry explicitly broken to U(l). Such a theory 
is described by the Lagrangian density: 

1 /If 2 m 2 h 

c = -d^d^ + — $t$ + — $t T3 $ _ _($t$)2 (20) 

where $ = (<&i, <&2) is a complex scalar doublet and r 3 is the 2x2 Pauli 
matrix. After rescaling as in equations (2)- (4) we obtain the equations of 
motion for $12 

fyd"$i,2 " (" 2 ± l)$i,2 + ($ f $)$i,2 = (21) 

where the + (-) corresponds to the field $1 ($2)- 
Consider now the ansatz 

with boundary conditions 

lim/(p) = 0, hmg'(p) = (23) 

\im f(p) = {a 2 + lf/\ }^9(P) = (24) 
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Vortex: with Non — Symmetric Core, a = 2.8 
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Figure 7: Field configuration for a non- symmetric- core global string with 
a = 2.8. 

This ansatz corresponds to a global vortex configuration with a core that 
can be either in the symmetric or in the non-symmetric phase of the theory. 
Whether the core will be symmetric or non-symmetric is determined by the 
dynamics of the field equations. As in the wall case we expect the existence 
of a critical value for the parameter a such that for a < a cr the vortex 
core is in the symmetric phase (g(p) = 0) while for a > a cr the symmetric 
core configuration is unstable towards decay to a new vortex configuration 
with non-symmetric core (g(0) ^ 0). To determine the value of a cr we 
solved numerically the system (21) of non-linear complex field equations with 
the ansatz (22) for various values of the parameter a. We used the same 
relaxation technique discussed in the previous section for the case of walls. 
For a < a cr ~ 2.7 the solution relaxed to a lowest energy configuration with 
g(p) = everywhere corresponding to a vortex with symmetric core (Fig. 6). 

For a > a cr ~ 2.7 the solution relaxed to a configuration with g(0) ^ 
indicating a vortex with non-symmetric core (Fig. 7). Both configurations 
are dynamically and topologically stable and consist additional paradigms of 
the defect classification discussed in the introduction. 

The existence of a defect with non-symmetric core corresponding to each 
defect with a symmetric core and winding at infinity can be established for 
various other cases. For example, our discussion for global vortices with non- 
symmetric core can be trivially generalized to the case of global monopoles. 
Another interesting extension is the case of non-symmetric gauged defects 
like Nielsen-Olesen vortices or gauge monopoles. Extension of our work in 
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these directions is currently in progress [|S . 

A common feature of the defects we have discussed so far is the fact that 
they emerge in theories where a larger global symmetry is partially explicitly 
broken. This is not always necessary. In the case of semilocal strings a U(l) 



gauge symmetry is embedded in a SU(2) global one|33|. The Lagrangian 
density is given by 

C = --F^F^ + 1 -{D^)\D^) - vm (25) 

with $ = ($ 1 ,$ 2 ) a complex doublet, D^ = <9 M — ieA^ a U(l) covariant 
derivative and V(<&) = ~($T$ — rf) 2 i.e. the 577(2) global symmetry is 
spontaneously broken with the U(l) part of it being gauged, the 577(2) 
symmetry is not explicitly broken, instead only the U(l) part of it is gauged. 
This is an alternative way to single out an S 1 subspace of the S 3 vacuum 
manifold. A similar analysis as the one discussed above (solution of a coupled 
system of non-linear ode's by a relaxation method) and an ansatz of the form 
(22) with A^ = eev(p) leads to a critical value for the parameter (3 = ^ (the 
only parameter of the model) such that for (3 < /3 cr = 1 the vortex with a 
symmetric core is stable. This result is in agreement with previous studies [J37], 
For f3 > (3 cr the symmetric core vortex is unstable to a configuration 



with a non-symmetric core called the 'skyrmion' |H], [3(], p9|1 . However in 
this case, since there is no explicit breaking of S77(2), the 'skyrmion' is not 
topologically stable (the vacuum is S 3 not S 1 ) and energetics favor expansion 
and eventual decay to the vacuum. 

4 Conclusion 

We have studied the existence and stability properties of topological defects 
with non-symmetric core and non-trivial winding at infinity. These defects 
possess cores which are either symmetric or non-symmetric depending on 
the range of parameters of the theory. They are distinct from previously 
discussed unstable hybrid f26|, |] defects as they are stable and are not nec- 



essarily connected to any other type of defect formed in a previous phase 
transition. 

These results have several interesting implications and potential exten- 
sions. For example, the observed non-spherical collapse of wall bubbles with 
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non-symmetric core may imply that the domain wall network simulations 
need to be re-examined for parameter ranges where a non-symmetric core in 
energetically favored. 

Our results are also relevant for baryogenesis mechanisms based on topo- 
logical defects|49|, ^(J, [51f| . These mechanisms are based on unsuppressed 
B+L violating sphaleron|52| transitions taking place in the symmetric core 



of the defects during scattering processes |53|, |54], [55], q6[ . Thus, the demand 
for symmetric core in such scenaria may lead to constraints in the parameters 
of the corresponding models. 

The existense of gauge strings and monopoles with non-symmetric core 
is also an interesting extension of our results. Do monopoles exist with an 
SU(3) x SU(2) x £7(1) symmetric core and are they more stable than the 
fully symmetric SU(5) one? The microphysical and cosmological implica- 
tions of defects with non-symmetric cores is an interesting area for further 
investigation. 
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